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SUMMARY 


«/-?*  K* 

^ Measurements  have-been  made  of  the  wall  pressure  fluctuation  and 
the  response  of  the  panel  structure.  The  plate  ^s,,  jMccited  by  a 
boundary  layer  (M=2L03)  and  shock  inpingement  angle)  the  f  low 

downstream  of  the  shock  is  separated.  The  wall  pressure  fluctuations 
and  panel  displacement  are  notably  greater  than  the  case  without  the 
shock  present.  ’ 

Theoresp6nse  of  a  plate  clamped  in  a  baffle  is  determined.  The 
plate  is  coupled  acoustically  to  a  uniformly  moving  fluid.  The 
boundary- layer  exciting  the  plate  does  not  interact  with  the  acoustic 
field  and  is  furthermore  taken  to  be  spatially  uncorrelated. 

Finally,  the  radiation  is  determined  over  the  solid  angle  of  the 
Mach  cone.  The  response  is  derived  using  finite  Fourier  transforms 
and  then  solving  the  resulting  algebraic  equation  by  expressing  the 
solution  expanded  in  a  complete  set  of  functions,  a  process  which 
leads  to  an  infinite  set  of  linear  equations.  The  right-hand  side 
of  this  system  of  equations  is  a  stochastic  variable  so  that  per¬ 
forming  an  ensemble  average  the  resulting  equation  is  solved  for 
the  covariance  matrix  which  is  subsequently  diagonalized  to  yield  the 
power  spectral  density  in  the  statistically  independent  states 
(degrees  of  freedom) .  ’  j  . 

Description  of  the  actual  computations  will  be  reported  at  a  later 
date. 
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1.  INTRODUCTION 

A  large  aircraft  la  supersonic  flight  undergoes  large  variations  is  flow  field 
over  its  surface.  This  paps?  is  concerned  with  studying  die  response  of 
a  structure  excited  fay  oonveeted  turbulence  at  nearly  aero  pressure  gradient 
and  by  shock-boundary  layer  interaction,  with  die  inclusion  of  the  coupling 
due  to  the  acoustic  field  on  each  side  of  a  panel.  Shock  waves  on  thin -walled 
structures  can  Impose  severe  loading  problems,  the  most  common  of 
which  is  the  self-induced  oscillation  which  is  generated  fay  an  oscillating 
shock.  The  shock  wave  can  easily  couple  with  the  forcing  frequency  present 
in  the  environment,  including  panel  resonances. 

From  interior  noise  point  of  view,  die  upper  region  of  the  airplane  fuse¬ 
lage  is  considered  the  principal  noise  radiator.  The  aerodynamics  in 
this  region  are  known  from  the  Prandtl-Mayer  relation,  and  further  down¬ 
stream  by  shock-boundary  interaction.  In  addition,  the  fuselage  skin  ex¬ 
periences  traveling  shock  waves  which  run  up  and  down  toe  skin  during  toe 
acceleration  period,  which  might  last  twenty  minutes  for  a  Mach  3  airplane. 

la  supersonic  flight,  the  vibration  of  the  surface  is  influenced  by  the  back 
pressure  resulting  from  the  radiation  of  sound  on  both  sides  of  toe  surface, 
so  that,  the  surface  motion  and  radiation  are  coupled  phenomena.  The 
interior  noise  level  is  determined  by  skin  panel  vibrations.  For  radiation 
below  toe  critical  frequency,  the  major  source  of  sound  arises  from  the 
interaction  of  the  bending  wave  with  toe  discontinuity  of  the  boundary. 

Above  the  critical  frequency,  toe  action  of  discontinuities  like  tear  stoppers, 
etc. ,  have  little  effect  on  altering  sound  pressure  level,  since  the  sound 
radiated  by  the  panel  is  in  the  form  of  Mach  wave  radiation. 

The  experiment  described  in  this  paper  indicates  that  some  simplifications 
in  the  model  can  be  made,  viz.  (1)  that  there  is  no  significant  interaction 
between  toe  plate  and  the  aerodynamic  forces  on  the  plate;  and  (2)  that  toe 
panel  displacement  is  small  in  comparison  to  its  thickness  so  that  thin 
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plat©  theory  may  be  used.  The  plate,  is  however,  acoustically  coupled  to 
the  external  flow  field  and  the  internal  cavity. 

Lyamshev  (1868)  has  solved  a  similar  problem  for  a  complex  structure. 
Dowell  (1869)  computed  the  transient,  non-linear  response  of  a  simply 
supported  plate  coupled  to  an  external  flow  field  and  a  cavity.  D&ygadlo 
(1867)  presented  a  linear  analysis  allowing  mutual  interaction  between  the 
plate  and  the  external  flow.  Fahy  and  Pretlove  (1987)  have  computed  a  first 
order  approximation  to  the  acoustic  coupling  of  a  flexible  duet  wall  to  the 
flow  field  through  the  duct.  Maidanek  (1886)  considers  an  infinite,  orthotropic 
plate  coupled  acoustically  to  an  external  flow  field.  Numerous  other  investi¬ 
gations  have  been  reported  on  acoustically  coupled  structures  with  varying 
degrees  of  approximation,  Irgens  and  Brand  (1868),  White  and  Cottis  (1968), 
Strawderman  (1967),  Creighton  (1970),  FfOwcs-Williams  (1966),  Crighton 
and  Ffowcs-Williams  (1969),  Dolgova  (1969),  Feit  (1966),  Lapin  (1967), 

Pal'tov  and  Pupyrev  (1967). 


MEASUREMENTS 

a)  The  Experimental  Arrangement 

The  flow  investigated  was  the  sidewall  boundary  layer  of  the  Jet  Propulsion 
Laboratory  20-inch  supersonic  wind  tunnel:  the  shock  was  induced  by  a 
30*  wedge  mounted  outside  the  boundary  layer,  off-center  and  on  the  same 
side  that  the  measurements  were  made.  This  was  done  to  offset  the  position 
of  the  reflected  shock  from  the  opposite  wall.  The  position  of  the  shock  was 
determined  by  observing  the  displacement  of  a  line  of  tufts,  and  by  a  static 
pressure  survey.  For  zero  pressure  gradient  detail  of  flow  field  and  panel 
response  has  been  previously  reported  by  Maestrello  (1968). 

The  experiment  was  arranged  to  perform  three  basic  measurements: 
mean  velocity  profile  ahead  of  the  shock  with  static  pressure  distribution 
across  the  shock,  wall  pressure  fluctuations  and  measurement  of  displace¬ 
ment  response  of  a  simple  panel  structure.  The  titanium*  test  panel 
measured  12  x  6  x  .  062  inches  and  was  brazed  on  ail  four  sides  of  a  3/4 
inch  x  3/4  inch  titanium  frame.  The  brazing  was  intended  to  simulate  the 
clamped  edge  condition.  The  panel  formed  most  of  one  wail  of  a  rigid 
cavity  measuring  14  x  8  x  6. 6  inches.  The  other  surface  of  the  panel  was 
exposed  to  the  flow  field.  The  pressure  differential  across  the  panel  was 
variable.  The  experiment  was  conducted  at  two  pressure  differentials, 
viz,  0.06  and  14  psi;  the  latter  corresponds  to  the  actual  differential 
between  wind  tunnel  pressure  and  local  ambient. 

The  side  wall  of  the  tunnel  was  modified  to  accommodate  two  identical, 
rigid,  steel  plates,  which  supported  the  necessary  instrumentation.  One 
plate  contained  an  array  of  holes  in  which  pressure  transducers  were 
mounted.  The  pressure  transducers  were  mounted  on  the  center-line 
of  the  tunnel  in  the  streamwise  direction  at  the  same  locations  where  the 


*TI-6AL-4V  Titanium  alloy  containing  6%  aluminum,  4%  vanadium,  90%  titanium 
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mean  static  pressure  measurements  were  made.  Two  types  of  pressure 
transducers  were  used;  one,  the  conventional  lead-zirconate  titaaate  type 
made  by  Atlantic  Research,  die  other  a  capacitance  type  made  by  Photocon 
Corporation  with  sensitive  diameters  of  0.06  inch  and  0.00  inch  respectively. 
Correction  due  to  finite  size  transducers  was  made  adopting  the  Corcos 
(1963)  approach.  The  panel  displacement  was  measured  with  Photocon 
capacitance,  displacement  transducers  mounted  on  brackets  which  could 
slide  along  a  bar  and  could  be  set  precisely  by  means  of  a  screw  mechanism. 

The  output  of  both  pressure  transducers  and  displacement  transducer  were 
recorded  oa  Ampex  FR-X800H  14-channel  tape,  recorded  in  the  FM  mode. 
Four  channels  were  used  for  simultaneously  recording  data  for  correlation 
measurements.  The  maximum  dynamic  range  was  obtained  by  splitting 
each  data  channel  into  two  tape  tracks  throu^i  phase  matched  filters  to 
separate  the  lower  and  higher  frequencies. 

b)  The  Wall  Pressure  Field 

Measurements  indicated  that  the  flow  field  in  front  of  the  shock  closely 

approximated  the  properties  of  equilibrium  of  an  adiabatic  flat-plate  boundary 

layer  (Maestrello  1988).  The  flow  in  front  of  the  shock  has  the  following 

characteristics:  Mach  number  M  =  3.03,  free  stream  velocity 

€ 

U  =2,100  ft/sec,  total  temperature  T,  =  567°  R,  boundary  layer  thickness 
©  * 

8=  1.37  inch,  boundary  layer  displacement  thickness  8*  -  0.445  inch, 

momentum  thickness  =  0.083  inch,  Reynolds  number  R  =  U  §/U  = 

5  ~3  -  - 

4.87x10  ,  skin  friction  coefficient  =  1,27  x  1C  ,  and  R^=  39.8 

Coles  parameter  (Coles  1964), 


The  pressure  ratio  across  the  shock  is  a  well  defined  function  of  Mach 
number,  for  a  15*  half-cone  angle,  the  pressure  ratio  is  approximately 
8.5.  Experimental  results  show,  however,  that  this  ratio  is  considerably 
smaller  (Ap  =  2.3).  It  is  postulated  that  interaction  with  an  expansion 
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wave  originating  at  the  base  of  file  wedge  Is  responsible  for  lowering  the 
pressure  differential  and  prcckclng  an  effective  decay  downstream, 

Figure  1c  In  the  present  case,  the  wedge  angle  induces  a  shock  in  the  boundary 
layer  large  enough  to  cause  a  separation:  farmer  downstream,  the  flow  be¬ 
comes  reattached  and  goes  back  to  the  flat  plate  condition.  This  transition 
takes  place  within  a  few  boundary  layer  thicknesses. 

Downstream  of  the  shock,  the  ratio  of  the  mean  pressure  distribution 
p  ,/p  and  the  ratio  of  the  rms  pressure  fluctuation  p'  ,/p'  vary  with 

Su  3  SO.  6 

a  consistent  relationship  and  both  reach  a  maximum  at  x/8  —  2.3,  where 
subscripts  s  and  sd,  mean  upstream  and  downstream  of  the  shock, 
respectively,  Figure  1.  Beyond  x/8  **  8  the  effect  of  the  shock  on  the  static 
pressure  vanishes.  Kistler  (1963)  indicates  a  similar  behavior  between 
mean  and  fluctuating  pressure  in  the  separated  region  ahead  of  a  forward¬ 
facing  step  at  the  same  Mach  number  and  upstream  Reynolds  number.  The 
differences  in  the  flow  geometry  only  alter  the  magnitude  of  the  pressure, 
in  that  the  ratio  of  the  mean  pressure  to  the  fluctuating  pressure  Pg(j/P’S(j 
14  in  the  present  experiment  while  Kistler  found  that  pg^/  p'^— 32„ 

The  normalized  power  spectral  density  measured  upstream  and  downstream 
of  the  shock  are  shown  in  Figure  2.  The  spectra  are  normalized  by  requiring 

V'CO 

/ ir  (03)  d  =  1  in  order  to  demonstrate  the  deviation  from  the  zero  pressure 

gradient  ease.  For  the  spectra  just  downstream  of  the  shock  more  energy 

is  concentrated  in  a  narrow  low  frequency  band  while  further  downstream  at 

x/6^4,  the  energy  is  distributed  over  a  much  broader  bandwidth  and 

approaches  the  shape  and  level  of  the  spectrum  taken  upstream  of  the 

shock.  The  normalized  power  spectral  density  found  upstream  of  the 

shock  corresponds  to  the  zero  pressure  gradient,  and  peaks  at 

<•>8  /U  —  2  while  downstream  the  spectral  density  is  modified  in 
© 

the  region  below  the  peak.  It  is  significant  that  by  altering  the  local  flow 


Figure  2.  Power  Spectral  Density  of  the  Wall  Pressure  Fluctuations 
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conditions,  only  the  low  frequency  ends  of  the  spectra  are  appreciably 
affected.  It  is  noticed  that  the  pressure  fluctuation  measurements  at  x/$ 

=*0  where  the  shock  impinges  show  a  noticeable  deviation  from  the  general 
pattern  in  the  higher  frequencies.  This  is  attributed  to  an  intermittent 
signal  superimposed  on  the  regular  pressure  signal  as  seen  on  the  oscillo¬ 
scope.  It  is  possibly  due  to  the  characteristic  fanning  of  the  shock  as  it 
goes  through  the  boundary  layer. 

Measurements  of  the  cross-correlation  are  shown  in  Figure  3.  The  cross- 
correlation  characteristics  are  a  function  of  position  downstream  of  the 
shock.  The  cross-correlation  between  positions  x/8  =  0.33  and  x/8  =  3.80, 
the  farthest  apart, has  characteristics  similar  to  those  found  at  zero  pres¬ 
sure  gradient  boundary  layer  in  that  the  ratio  between  the  convection  velocity 

and  the  freestream  velocity  U  /U  =  0.72  and  that  the  correlation  between 

C  6 

those  two  points  is  still  significant.  The  cross  correlation  of  the  shortest 

distance  between  x/8  =  0.33  and  x/8  -  0.75,  shows  that  the  convection 

velocity  is  very  low  U  /U  =0. 13  and  the  correlation  i6  very  weak.  The 

C  6 

correlation  between  x/8  *  0.33  and  x/o  =2,25,  where  x/8  =2.25  cor¬ 
responds  to  the  maximum  static  pressure  ratio  is  negative.  The  shock 
induces  the  boundary  layer  to  separate  and  the  recirculation  within  the  separa¬ 
tion  region  permits  the  sign  of  the  pressure  to  change.  Kistler  argued  that 
the  fluctuating  pressure  in  the  separated  region  arises  from  the  combined 
action  of  the  turbulent  shear  layer  and  the  recirculating  flow.  The  picture, 
however,  is  not  yet  clear  enough  to  develop  a  model  for  time  dependent 
loading,  since  the  geometry  of  the  separated  region  is  the  primary  variable 
in  estimating  the  pressure  amplitude  and  resulting  phase. 

No  measurement  of  the  lateral  cross-correlation  was  made  during  the  test; 
however,  for  the  purpose  of  computing  the  response  of  the  panel,  it  is  assumed 
that  the  pressure  decays  similarly  to  that  in  the  case  of  zero  pressure  gradient 
e~  ^  *2  where  a  =  0.26  and  V  is  the  spatial  separation  (Maestrello  1968). 

Ct 
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This  choice  overestimates  the  lateral  cross-correlation,  since  the  flow 
field  is  far  from  being  homogeneous.  However,  the  overestimation  may 
not  be  exceeded  by  a  factor  of  2. 

c)  The  Panel  Response  Field 

Measurements  were  made  of  the  power  spectral  density  and  cross-correlation 
of  die  displacement.  Typical  results  are  shown  in  Figures  4,  and  5  for  a 
pressure  differential  of  14  psi.  The  static  deflection  of  the  panel  was  0. 06 
inches  at  the  center,  and  the  dynamic  deflection  was  small  in  comparison 
with  its  thickness. 

The  displacement  spectral  density  at  the  center  of  the  panel  show  pronounced 
spikes,  the  lowest  frequency  of  which  corresponds  to  the  lowest  mode  of  the 
panel.  The  accuracy  beyond  a  frequency  of  3100  Hz  was  poor  due  to  the 
spatial  resolution  of  the  capacitance  transducer,  and  therefore  the  spectrum 
beyond  3100  Hz  was  ignored. 

Space-time  correlation  measurements  were  made  along  the  panel  centerline 

from  x  =  x'  =  3  in.  y  =  y*  =  3  in.  at  one-inch  intervals  up  to  a  maximum 

separation  of  6  in.  The  correlogram  indicates  a  convected  feature  with 

a  phase  velocity  ±  U  =  770  ft/sec.  This  convection  velocity  corres- 
cp 

ponds  to  that  found  in  the  previous  experiment  using  the  same  arrangements, 
except  that  no  shock  was  present  (Maestrello  1968). 

In  comparing  the  results  of  the  present  and  previous  experiments,  it  is 
concluded  that  the  sign  change  of  the  convection  velocity  is  attributed  to 
the  presence  of  the  shock.  Furthermore,  the  cross-correlation  of  the 
wall  pressure  also  reflects  a  phase  change  for  a  separation  of  2.5  inches, 
which  is  in  the  same  location  as  the  phase  change  which  occurs  for  the 
displacement  correlation  in  Figure  5. 
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Figure  4.  Displacement  Spectral  Density 
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Figure  5.  Broad  Band  Space  Time  Correlation  of  the  Panel 

Displacement  Along  the  Center  from  x  =  xT  =  0.25  Ft. 
y  =  y*  =  0.26  Ft. 
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S.  ANALYSIS  OF  ACOUSTICALLY  COUPLED  PANELS 
a.  Two-dimensional  Finite  Panel 

The  vibration  of  the  panel  is  induced  by  an  arbitrary,  external  pressure 
field  F.  It  is  assumed  that  the  panel  motion  does  not  interact  with  tke 
turbulent  boundary  layer,  i.  e. ,  the  forcing  field  is  not  altered  by  the 
plate  motion.  However,  the  panel  is  acoustically  coupled  to  the  fluid  on 
both  sides  of  the  panel, 

The  equation  of  motion  for  an  harmonic  component  of  the  displacement, 

W,  of  a  thin  panel  with  a  force,  F,  and  a  pressure  differential,  P2  “  Pj  +  <$P 
acting  upon  it,  obeys  the  equation 

B  A2  W  -  pp  u  W  =  F  +  p2  -  p1  +  Sp  (1) 

where  the  bending  stiffness,  B,  may  include  hysteretic  damping,  and  where 
p  is  the  mass  per  unit  area  of  the  panel, « is  the  angular  frequency,  p0 
is  the  acoustic  pressure  on  the  streamside  of  the  panel,  p.^  is  the  acoustic 
pressure  below  the  panel  and  g  is  the  static  pressure  differential. 

The  perturbation  pressures,  p  and  p  ,  are  related  to  the  velocity  poten- 

1  « 

tials,  which  satisfy  time-independent  wave  equations  in  the  appropriate 
regions.  In  solving  these  equations  one  uses  a  boundary  condition  which 
relates  the  potentials  to  the  panel  displacement.  These  relationships  may 
be  made  more  obvious  through  the  use  of  Green's  theorem.  Thus,  it  is 
required  to  solve  a  system  of  three  coupled  partial  differential  equations, 
the  first  of  which  is  not  separable  for  the  clamped  edge  boundary  condition. 

Pj  and  p2  may  be  found  directly  as  function  of  W.  Thus,  consider  first 
the  cavity.  The  acoustic  velocity  potential,  9,  satisfies  the  Helmholtz 
equation 

A  9  +  k2  9  =  0  (2) 


t 
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with  boundary  condition 
am 

where  =  ~i  o>W, 

Hie  Greens  function,  g,  for  a  cavity  with  hard  walls  satisfies  the 
equation 


-|2=  0  on  all  walls  except  on  the  plate 


A  g  fr  i  5*)  +k2g  F|r*)  =  -4tfS(r-r*) 
c 


(3) 


and  is  given  by  Morse  and  Feshbach,  Vol.  H  (1953) 


g  (r/r»)  = 


m7 rx  nwrx  n77y  nTry’ 

cos - cos - COS  COS  -r — — • 

A  AT  V  ^  &  w  w 

J  2  c _ c  c _ c_ 

a  b_  m=0  ufO  fcmfcn  k  sin  (k  d) 


c  c 


mn  mn 


(4) 


!cos 

cos 


where  k 


2 


mn 


k  z  cos  k  (zf  +  d) 
mn  mn 

k  z'  cos  k  (z  +  d) 
mn  mn 


-(f)'  M 


z  >  z’ 
z  <  z' 


and  k  =  <y/c  where  c  is  the  speed  of  sound  in  the  cavity  of  dimen- 
c  c  c 

sions  a  ,  b  ,  d* 
c*  c 


By  applying  Green’s  theorem,  the  integral  equation  for  9  is  obtained, 


_  i  [(  __  aw®  _ 
<r>  =  —  JJ  «*/?)  05~drs 


<P( r) 


4  7T 


(5a) 


plate 


Now  using  the  boundary  conditions,  this  becomes 


W  =  W  <?,d? 

plate 


(5b) 


The  pressure  is  related  to  by 


p1=-i<»po<p 

where  pc  is  mass  density  of  the  fluid  in  the  cavity. 

To  compute  p  ,  it  will  be  more  convenient  to  operate  with  the  differential 

u 

equation.  Let  the  acoustic  velocity  potential  in  the  flow  field  be  denoted 
by  ’A.  By  applying  the  Fourier  transform  on  the  (x,  y)  coordinates, 
one  gets  the  ordinary  differential  equation 


where 


+  C2  $  (<*»($,  z)  =  0 


C2  =  k2  +  (M2  -  l)a2  -  2  k  M«  -jg2 


(6) 


\fr(x,  y,  z)  =>J  j  dad/?ei(ax  z) 

-00  -CO 

K  =  w/c,  M  is  the  flow  Mach  number  and  c  the  speed  of  sound  in  the 
region  above  the  plate.  Only  the  positive  exponential  solution  to  equation 
(6)  is  chosen,  since  it  is  the  solution  representing  outgoing  waves. 

Thus, 

J  (a,0,  z)  =  A  (a,/3)  e^Z  (7a) 

The  boundary  condition,  arising  from  the  continuity  of  normal  displacement 
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is 


A 


UMlH  ,  ~ie  & 
da 


where  the  differential  operator 


L  -  ktiMf 
ox 


Thus, 


A  A  . 

^  ,  o  ,  LW  i£s 
y  (<*, 8,  ^  =  -cj8 


(7b) 


Now  since 


LW  =J  f  dx’tfy1  e“l^X*  "^y\w  0c1,  yf) 


o  o 


then 


a  b 


»//  (x,  y,  a)  =  -  -\jj G  y.  z  I  x‘»  y*»  °)LW  (x!,  y!) 


@) 


0  o 


where 


»  «>  i  [a  (x  -  x’)  +/?(y  -  yJ)  +  <(  z  -  z‘) 

G  (r  |  ?)  ~  j  j  dad/3  5 j -  (9) 


-CO  -00 


which  is  found  in  Appendix  A  to  be  for  supersonic  flow, 


2iri  e 

M7Ti 


G  (r  |  r!) 


i  K  (Mu  +  V  u2  ~ 


vrrz 


o 


outside  the  Mach  cone 
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and  for  subsonic  flew. 


G  (r  1  r’>  = 


2m  e 


iK(Mu  +VT3, 


Vl  -  M2  Vp**H2 


X  -  X1 


except  In  this  case,  k  =  V^_  m2  ^  u  =V71 


A 

If  LW  had  been  evaluated  as 


A 

fle  -ft  M)  W 


(10) 


then  equation  (8)  would  read 

a  b 


y,  z)  =  J J  dx’dy*  W  (xJ,  y’)  I?G  (x,  y,  z  jx?,  y,  z *)  (11) 

47r  o  o 

This  equation  is  formally  correct  if  L*G  Is  interpreted  as  a  distribution, 
which  is  to  say  that  one  partially  integrates  to  obtain  equation  (8)„ 


Now  using  equation  (8) 


P2  (x,  y,  z)  *  -ip0c  L^(x,  y,  z) 


lP0° 


2  a  b 


4  ft  o  o 


f  f  dx’dy1  G  (x,  y,  z  jx1,  y\  o)|lF  W  (xs,  ys) 

(12) 


where  po  is  the  density  of  the  fluid  above  the  plate,  and  where  partial 
integration  has  been  utilized.  Had  equation  (11)  been  used  instead, 
equation  (12)  would  read 

ip  c2  a  b 

P2  (x,  y,  s)  =  J J' dx'dy5  W  (x,  y)  Ilf  G  (x,  y,  z  j  x1,  y},  o)  (13) 


O  0 


5 


which  is  reducible  to  equation  (12)  by  partial  integration.  Thus,  supersonic 
flow  does  not  present  any  especial  difficulty  aside  from  the  fact  that  G 
is  singular  all  along  the  Mach  cone,  and  this  is  an  integrable  singularity. 

Inserting  the  expressions  for  pg  and  into  equation  (1)  results  in  a 
single  partial  integro-differential  equation  to  solve,  viz. , 


ip  c2  a  b 

BA2W-  P  <u2w=F  +  S  +— —  f  f G(x  ,  y,  |  x>,  y'jllf  W  (x!,y!)dx*dyl 
P  P  47 T  iJ0 


■a  Pc  //  g  (xc,  yc  |X^,  y*)  W  (X,  y)  dx  dy 

4  7T  v  ^ 

plate 


(14) 


where  the  subscript  c  refers  to  the  cavity,  thus 

a  -  a 
c 


Equation  (14)  presents  a  formidable  computational  problem.  The 
Green's  function  g  is  known  as  an  infinite  series  which  is  slow  to 
converge  (l/n)  thus  compounding  the  difficulty  by  an  increasing 
number  of  necessary  operations  to  maintain  a  given  accuracy. 

An  alternative  to  solving  equation  (14)  is  to  convert  it  to  an  integral 
equation  for  its  Fourier  amplitudes  and  to  solve  the  resulting  equa¬ 
tion.  The  advantage  is  that  this  equation  is  simpler  (though  it  is  a 
singular  integral  equation).  The  following  notation  shall  be  employed: 


A  ->.  1 

f  (K)  =  -r- 

2  Tt 


supp  Jf| 


,  -*  “iK  •  r  ,  .•*•, 
are  f  (r) 


f<?)'2?-  ^  dS  e‘K  rf® 
K- space 


The  result  of  applying  the  Fourier  transform  to  equation  (1)  is 


?'  2  A  A  A  A 

BA  w-  p  0)  W  =  f  +  p„  -  p 
p  -  1 


where 


f  =  F  +  5 


The  first  term  in  (15)  may  be  evaluated  using  Green's  theorem;  thus, 
A2  W  =  K4  W  +£d  ?g  je_i^  °  r  ( A  W  +  K2  W)] 


Now,  for  a  plate  clamped  on  its  edge  to  a  rigid,  plane  support,  the 
following  boundary  conditions  hold, 


a2  W 


d  W 
n 


a3  w 
an  a  s2 


on  edge 


Where  s  is  in  the  direction  of  the  edge,i.  e. ,  the  tangent,  Thus, 


A2  W  =  K4  W  +  A  [w] 


(16b) 


where 


If  more  general  boundary  condition*  are  to  be  oonsidered  (e.g,  elastic 
foundation)  the  above  expression  must  be  replaced  by  the  right  side  of  18a. 
From  the  equation  (7b)  and  the  equation  prior  to  equation  (12)  It  is 
found  that 


£  a  lP0°2  <k-«M)2  & 

2  ■' 


a  * 

w& 


(17) 


From  equation  (5b)  it  is  found  that 


v  i?.»w ff-> 


slate 


=  -0[w] 


where 


g  <£l 


l  #  - 1  k  •  7  ,  4., 

r)  J  dr  «  g  (r|rf 


plate 


(18) 


Substituting  these  results  Into  equation  (15)  gives 


T& w£)  + 


BA 


W 


(19) 


where 
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Let  denote  the  finite  Fourier  transform  of  a  beam  eigen  function, 
9^.  It  follows  from  the  Fourier  representation  that  the  form  an 
orthogonal  set  on  the  infinite  interval.  Thus,  expanding  W  as 

=  2  Wmn  <“  “>  ^b> 
m,  n 

or  alternatively,  W  as 

W  (r)  =  y  W  <p  (-}  (p  u) 

^  mn  Tm  \a/  n\b/ 
m,  n  '  •  '  • 

and  introducing  these  expressions  into  (19)  and  subsequently  utilizing 
the  orthogonality,  gives 


(20) 


(21) 


w  +  2  r  w  =  & 

mn  r,  s  mnrs  rs  ^mn 


(22) 


where  <J)  are  the  projections  of 
mn 


A  -o 


*2  4>mn 
m,  n 


and 


(See  Appendices  B ,  C  and  D) 


! 
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The  computation  of  the  integral  F  may  be  simplified  by  deforming  the 

mnrs 

contour  on  the  oc-  plane.  IXie  to  the  manner  in  which  the  Fourier  transform 
was  chosen,  the  integrand,  except  the  term  T(ic),  is  single-valued  and 
analytic  in  fee  lower  half -plane.  The  contour  will,  thus,  be  deformed  in 
this  half-plane.  This  deformation  is  determined  by  the  analytic  properties 
of  the  function  T(k),  equation  19  a.  2 

2  /o«  /,  a2M2\ 

a.fl)  -  («2  +  j82)2  - i_A - JsL 1 


T(a.]8) 


B 


\/  ,  2  ~2  ~  2  ^2 
v  k  +  (M  -l)a  -  2  k  Ma  -  /? 


This  function  is  two -sheeted  wife  square-root  type  branch  points  at 


kM  *V  k2  +  (M2  -l)#2 


M2  -  1 


The  sheet  associated  with  fee  positive  value  of  fee  squre  root  will  be 
termed  the  physical  sheet,  since  it  corresponds  to  outgoing  radiation. 

The  function  has  ten  zeros  on  the  two  sheets,  four  zeros  on  each  sheet 
with  the  same  values,  corresponding  to  resonances  of  the  plate  and  the  other 
two  zeros  are  located  near  the  branch  points  on  one  of  fee  two  sheets,  inde¬ 
pendent  of  each  other. 


It  is  convenient  to  make  the  following  substitutions 


M  =■ 


4  Pn“> 

and  Y  =- 


B 


The  equation  (19a)  may  be  written 

T(a  >($)  =  («2  +j32)2  -  Y4 


Vk2  +  (m2  -  1)  «  -  2  k  M«  -  j32 


nsfr  jt  .-*t^*- 
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In  the  present  case  fi  is  a  small  number  (»  .  0015)  so  that  approximate 
values  of  the  zeros  may  be  found,  expressed  as  a  power  series  in  fi. 
To  the  second  order  th§pe  zeros  are 

^M2  2’ 


a*  =  A*  -it 
8  8  2 


m 


V  +(M2-1)  /32)(Af  +  2 If  02  +  j34  -  y4)2 


where 


A 


=  k  M  iV  k2  +  (M2  -  1)  fl2 


M2  -  1 


a 


=  A 


(±)l(i)2  (±)1(±)2  4  A 


mA,  .  ,  „  \2 

iY4 

\  *  1 

i 

(i,X<*>2  <PH)l(i)2>  IV  +  <m2-1)  klhl%  2  k 


where 

aw1w2  *  t^V-/}2  (±)2  Y2 

The  last  four  zeros  exist  on  both  sheets.  The  location  of  the  first  two  zeros 

may  be  distinguished  into  three  possibilities:  when  Y<  A  then  a 
+  s  s 

and  a  are  respectively  on  the  unphysical  and  physical  sheets ;  as  Y  is 

6  + 

increased  such  that  A  <  Y  <  A  ,  then  a  moves  off  the  unphysical  sheet 

S  S  8  ^ 

and  crosses  over  to  the  physical  sheet  and  a  remains  unchanged;  as  Y 

is  further  increased  such  that  A  <  Y  then  a  crosses  over  to  the  unphysical 

s  s 

sheet  and  a  remains  unchanged.  A  typical  configuration  for  the  poles  and 
s 

branch  points  is  shown  in  Figure  6. 


-v— 


25 


The  function  T  has  been  analytically  continued  Into  the  k-plnne  by  giving 
k  a  small  negative  imaginary  part;  so  that,  the  branch  points  and  the  poles 
are  displaced  off  the  real  axis  as  indicated  in  the  previous  figures. 


The  branch-cut  integral  is  given  by 


00 


I(j3)  =  /dt 


/■ 


4>m  (a «)  [BA(qtj8)  -  Q  (a.fl)] 


W  +  /j 


pno2(k-  M)2 


2  2 
Pn  c  (k-«M) 


(M2-!^2  -2kMa-g2, 


(B  («2  +f?)2  -  P  «V  -Ury  - j - 

P  \y  K  +  (M  -  1)<*  -2k 


_y 

Ma-02  / 


where  a-  A  -it 
s 


So  that 


r. 


mnrs 


/*  /•» 

2w  i  2  residues  d  fi  +  J  (fib)  l{{3)  d(3 

-00 


The  branch-cut  integral  is  exponentially  damped  rather  oscillatory,  so  it 
may  be  readily  performed  numerically  using  Laguerre-Gauss  quadrature. 

2  7T 

The  second  integral  is  more  difficult,  it  oscillates  with  a  period  . 


In  solving  equation  (22)  maximum  values  of  the  indices  are  postulated  . 
This  is  justified,  since  the  index  is  inversely  proportional  to  some  length 
on  the  panel.  Now  there  certainly  exists,  from  the  experimental 
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point  of  view,  a  smallest  length  to  which  a  disturbance  may  be  localized. 
After  having  solved  equation  (22)  the  plate  displacement  is  simply  the 
Fourier  transform  of  (21),  thus, 

wff>= ^  UH(fr)  w 


where  <jp  is  a  beam  eigen  function, 
m 


To  find  die  sound  pressure  level  in  the  cavity  the  expression  for  W 
from  equation  (23)  is  inserted  into  equation  (5b)  to  give 


2 

r.  0)  P 

Px  (?  “  c 


mTtx 


nTiy^ 


a  b 
c  c 


2  2  w  i  j 

m,n  r,s  rs  mr  ns 


cos - *■  cos  ,  cos  k  (z+d) 

d  mu 


a 


_S_ 


sin  k 


where  I  and  J  are  given  in  the  appendix  B. 
mr  ns 

Similarly,  the  radiation  may  be  computed  from  equation  (12) . 

The  force  is  not  a  deterministic  function  as  has  been  implicitly 
assumed  from  the  outset,  but  a  stochastic  variable  whose  correlation 
properties  are  known,  either  via  a  model  or  directly  from  experi¬ 
mental  data.  Thus,  it  would  only  be  meaningful  to  compute 
statistical  averages  of  the  response  based  on  statistical  averages 
of  the  force  (i.  e.  cross-correlation) .  The  procedure  to  be  used 
is  an  amendment  of  a  procedure  due  to  Rosenblatt  (1962),  the  no¬ 
tation  is  that  of  Rosenblatt. 

Consider  a  homogeneous,  stationary,  random  process  X* 

r ,  t 

(Rosenblatt  writes  this  as  X*  ,(co)  to  explicitly  indicate  that  it 

r,  t 
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is  a  function  defined  on  &  sample  space)  with  a  cross -correlation  defined 
as 


B  (?,  t;  f) 


(24) 


where  <  )>  is  the  expectation  operator,  i.e. ,  R  is  defined  through  an 
ensemble  average.  Because  the  process  is  homogeneous  and  stationary, 


R  (r,  t;  ?',  f)  =  R  (r  -  i\  t  -  t’) 


or 


R(?-r»,  t-t*)drdt  =  <d  M  (x~  ^  d  M  (X-,  fc,)> 

where  d  M  (x  ♦  ^  is  the  Stieljes  measure  of  the  process.  The  procedure 

may  be  simply  stated  as  the  problem  of  finding  a  Fredholm  expansion  of 

R  and  subsequently  representing  X  by  such  an  expansion.  Such  an  ex- 

r,  t 

pension  is  provided  by  the  eigenfunctions  and  eigenvalues  of  the  integral 
equation 

03 

^(r,  t)  =  A  f  R  (?  -  r\  t  -  t*)  ^(r1,  t’)  dr'  dt 


(25) 


The  spectrum,  is  of  course,  continuous.  The  eigenfunctions  are  plane  waves 
and  the  eigenvalues  the  inverse  of  the  power  spectral  density  as  can  be 
seen  by  applying  the  Fourier  transform.  Thus  the  desired  expansion  for 
R  is 


R  (r  -  r»,  t  -t’) 


^  —  03 


Now  let 


Z-, 

K’<*> 


(2-77)  \/R  (£ 


=P 


(K  •  r  -«t’ 


Xr,  t> 


(26) 
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II  follows  from  (24)  and  (26)  that 

<zsuzeC«>  ’  «* -*>«<—'> 

Sine*  <6  M  {X-  jawfe  ,,)  >-  R  (?-?*.  t  -  f)  df  dt 

*t  •  *  *  * 

Thu*  ,  the  Z  arc  independent  random  variable*  with  unit  variance  and 

K, 

with  aero  mean  if  E  X*  »  0.  The  transform  of  (26)  is, 

ft  * 

xf  ,  -  /  VS  £,  u  >  ""’*!«  <»T) 

—  • 

A  simple  calculation  reveals  that  (27)  satisfies  (24), 

These  results  will  now  be  applied  to  the  plate  displacement.  Thus,  the 
cross-power  spectre!  density  (CPSD)  of  the  plate  response  is  given  by 
(asterisk  denotes  complex  conjugate) 

<W(?)WM?’£>-m  E  j  <Pm  (t)  %{i)  <WmaWr.”> 

(2«) 


Now  if  the  solution  to  equation  (22)  is  represented  as 


W 


mn 


Z 

i.) 


Yamij  ^ij 


then 


<W  W*  >  -  V 
mn  rs  '  Li 

ijkl 


mnij 


Yrskl 


<  $ 


*ki  > 


Further. 


w>=;jdKdKi 

-00-03 


_  &<«a)^{0b)tj'{fa'a)4/f  <0‘b) 
fJdKdK’~ - ‘ - T(ft)  f*  (K1) -  <f  (K)  f*  (K*>> 


The  force  F  is  now  identified  with  X-  A  so  that 

r,  t 


<f  {K)  f*(K*)>[$  (K,  w)  R  *  (K’,  W)]  <Zg  ^  Z|,^ 


[A  -  A  -»  11/2  . 

=  [R  (K,w)  R*  (K*,  W)  J  5(K-K') 


In  summary  then 


rm  _  (aa)  xf,  (jib)  |r{K,w)|»1/*  (aa) ip*.  (Rb) 

<w  w*  >=  2  y  4.  y*tl  I  dK- - 2 - - — - L  — 

nra  rs  mni>  ™  J  T(K)  2 


Analogous  to  (28) ,  the  expression  for  the  CPSD  of  can  be  written  in 

terms  of  E  W  W*  from  (29) .  The  same  can  also  be  done  for  the 
mn  rs 

radiation. 
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B.  One  Dimensional  Model 

To  simplify  the  computations  we  assume  that  the  transverse  plate  dimen¬ 
sion  is  very  large  and  that  no  flexural  waves  propagate  in  the  transverse 
direction.  With  these  simplifications  equation  (22)  may  be  written 


W  •  2  F  W  =  d> 

n  m  nm  m  n 


(30) 


where 


r 

nm 


/ 


-oc 


00 


dK 


^(Ka>  A[ym(l)] 


T  (K) 


where 


T(K) 


^k2  +  (M2  -1)  K2  -  2  k  K  M 


(31) 


and 


$ 

n 


(a  K)  F  (K) 
T(K) 


(32) 


The  circumvention  of  the  branch  points  in  the  above  expressions  will  be 
described  shortly. 


If  G  denotes  the  inverse  matrix  to  5  +  T  then  the  solution  to  (30)  is 

mn  ran  mn  ' 


W  =  2  G  <p 
n  m  nm  m 


Now  performing  the  ensemble  averages  as  before  gives 


31 


/W  W*  \  *  2  g  <f>*  \  G* 

N  m  n  /  r,  s  mr  \Tr  /  rs 


but 


/.«  \f/*  (a  K)  /*«  ik(aK‘) 

<+r ** >  - j d K~r^rr K' “F^  <F(K)F*(K'» 


To  make  the  discussion  concrete  let 


<F  (K)  F*  (K*)>  =  P  («)  6  (K»  -  (K  ) ) 

c 

which  corresponds  to  a  spatially  uncorrelated  pressure  field  with  convection 

velocity  U  and  power  spectrum  P  (to), 
c 


Thus, 


aw 


<^r^g  > 


P(<o)  /  dK 


•  K  (aK)^(aK-  -~) 


-03 


T(K)  T*(K-  £) 
c 


(34) 


The  major  contribution  to  the  integral  for  1^,  equation  (31),  comes 

about  when  the  peak  of  ip  is  close  to  the  peaks  of  1/  T  (K) ;  since  ip  is 

D  2*77*  v  ^ 

a  highly  oscillatory  function  (period  = — )  with  a  peak  at— —  and 

a  1  a 

decaying  with  the  distance  from  this  point  anc*‘Y(K)  is  a  non-oscillatory 

function  with  peaks  whenever  K  equals  the  real  part  of  the  pole3  which  are 

roughly  located  at  Y times  the  four  roots  of  unit  and  the  trapped  wave  poles 

k  k 

near  the  branch  points  - — -r— r  .  But  for  the  frequencies  we  are  con- 

1  +M  M  -1  ^ 

sidering  (up  to  3000  Hz)  only  the  pole  near  lies  on  the  physical 

sheet.  For  this  frequency  range  the  trapped  wave  pole  is  bounded  by  0  and  .  15 

v  m 

and  the  pole  near  Yby  0  and  .6.  Now,  the  peak  of  ip*  is  given  by - 

Q  ct 

which  is  numerically  (see  Appendix  B)  .155,  .26,  .36,  .465,  .57,  .67,  . 

,  -  13 

and  the  period  is  .206.  The  height  of  the  peaks  decays  roughly  asf 

is  .  05  of  the  first.  The  function  also  has  peaks  that  decay  in  the  same 

I  (K) 


'fe! 
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manner.  Thus,  the  infinite  matrix  has  appreciably  non-zero  entries 

only  in  the  upper  left-hand  corner .  Consequently,  we  need  only  compute 

T  for  the  first  4  or  -  modes,  say,  and  then  invert  this  matrix  +  I  to 
tan 

obtain  the  upper  left  square  matrix  of  order  4  or  8  of  G  and  thus  for 

ran 

higher  modes 


G 

mn 


So  that 


<WmW;>  - 

for  other  than  the  first  few  modes. 


The  contour  for  the  integral  in  (34) ,  Figure  8,  is  similar  to  one  described 

b 

earlier  but  the  term  T*  (K  -~)  introduces  further  poles  and  branch  cute. 
The  contour  and  ite  deformation  are  shown  in  the  Figures  8  and  9.  In 
Figure  9  only  the  contours  for  that  part  of  the  integrand  which  is  analytic 
in  toe  upper  half-plane  are  shown. 


\f/  (K)  may  be  written 


vL  (K)  =  S  (K)  +  T  (K)  e 
Ym  '  m  m 


where 


S 

m 


(K) 


i  a  (a  K)3 
m _ 

(aK)4- 


and 


Figure  9.  Deformation  of  contour  in  Figure  8. 
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~T  (K) 
2  m  '  7 


^[(Xm  +  1“maK(- 


COS 


chx 


m 


2  2 
» -  x* 


(a  K)2  + 


) 


-  (a  X  -iaK) 
'  m  Am  ' 


£ 


sin 


2  2 
’aK)  -X 
'  m 


(aK) 


77$ 

'  in 


with  and  X  defined  in  Appendix  B. 


Now 
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The  above  expression  is  evaluated  as  a  sum  of  residue  contributions  plus 
branch  -  cut  integrals. 
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Figure  10.  Location  of  poles  and  branch  cuts 


As  indicated  previously,  zg  and  zg  are  not  on  the  physical  sheet  for  the  range 
of  frequencies  considered  (Up  to  14000  Hz). 


where  T+  +iyj  is  the  expression  following  equation  (31)  the  plus 

being  included  only  to  explicate  that  it  is  the  positive  side  of  the  square  root 
branch  cut.  X  fcj  +-  ^  +  i  yj  is  obtained  from  T+  by  replacing  fi  by  -  fj. . 
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It  is  felt,  that  the  effect  of  the  pole  at  z  cancels  the  contribution  of  I  and 

O  1 

similarity  the  residue  at  z,_  cancels  the  contribution  of  I_.  Thus  only  I0 
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and  I,  will  contribute.  Therefore 
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By  the  same  argument  which  eliminated  I.  and  1^  we  may  replace  1%  and 

I  by  pretending  that  z  and  z  are  on  the  physical  sheet  and  thus 
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To  evaluate  the  residues  it  is  necessary  to  determine 
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In  a  similar,  though  simpler  fashion,  we  may  evaluate  the  integral  in 
equation  (32). 
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Thus  A  may  be  decomposed  into  factors  analytic  in  upper  and 
lower  half -planes,  respectively.  Denoting  these  terms  by  G  +  and  G~ 
we  have, 
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Figure  11.  Contour  for  F  equation  (SI) 


where 

G+  (K)  =  Am  (K)  ^Sn*  (K)  *  Tn*  (K)  e‘  K  *) 
and 

G“  (K)  =  Bm  (K)  (K)  +  Sn*  (K)  e“  1  K  a 

The  contour  of  integration  is  indicated  above  . 


The  contour  is  deformed  in  the  appropriate  haif-planes  and  branch  cuts  are 
replaced  by  poles  ad  discussed  previously.  Thus, 
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We  now  invert  the  matrix  S  +  T  .  Let  us  denote  this  inverse  by  Gti . 

mu  xnn  13 

Then 


/W  W  *\  =  G  <F  F  *\  G 
\  m  n  /  mp  v  p  g  /  gn 

where  the  tilde  denotes  Hermition  conjugate  and  sum  of  p  and  q  is  implied. 
The  final  step,  then,  becomes  the  diagonalization  of  the  CPSD  matrix  <^W  W  *} 
thus  giving  the  PSD  for  the  actual  degrees  of  freedom  of  the  system. 


C.  Acoustic  Power,  Power  Radiated 

In  computing  the  power  radiated,  we  make  use  of  the  unitary  matrix  which 
diagonalized  the  CPSD  matrix.  The  average  power  radiated,  to  the  far  field,  P, 
is  given  by 
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u  =  ±  j  (r-6) 

r  dr 


1 


42 


ik (vT 


M2  sin2  #  +  M  cos  # 


) 


M2  -  1 


The  pressure,  p,  a  micro  phene  in  the  far-field  would  measure  is  given  by 
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The  average  radiated  intensity  is  giver  by  the  expectation  value 
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and  where 
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Let  U..  denote  the  unitary  matrix  which  diagonalizes  <W  W  *\  ,  i.e. , 
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the  transformation  to  the  actual  degrees  of  freedom.  Further,  let  A 
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denote  the  power  in  the  m  degree  of  freedom,  then  we  have 
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The  integrals  have  already  been  encountered,  they  are  simply  the  finite  Fourier 
transform  of  the  9^  (x)  so  that, 
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APPENDIX  A 


If  we  make  the  following  changes  in  variables 
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then  expression  for  the  Green's  function  in  equation  (9)  becomes 
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The  contour  of  integration  for  the  above  integral  is  shown  below 

£  -  Plane 

Branch  Cut 
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I(u)  -  j  d  £  e1  ^ U 

By  considering  the  asymptotic  form  of  the  Hankel  function,  the  shore 
Integral  is  seen  to  be  convergent  in  the  upper  half-plane  for  values  of 
u  and  R  such  that 

u  >R 

that  is,  the  region  inside  the  Mach  cone 


The  contour  may  be  deformed  to  be  a  contour  along  the  branch  cut  as 
shown  below 
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We  thus  get  (see  Magnus  &  Oberhettlnger  p.  179) 
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COS  K  ~  Ch  K 
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and  the  normalization  is  given  by  (Dzygadlo  1967) 
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and  the  eigenvalues  are  the  roots  of  the  equation 
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and  are  approximately  given  by 
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J  is  the  same  as  I  with  a  and  a  replaced  by  b  and  b  ,  respectively, 
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APPENDIX  C 
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The  integral  may  be  found  by  appropriately  combining  the  following 
four  integrals 
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APPENDIX  D 
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and  is  given  in  appendix  B. 
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